Verification of a non-hydrostatic dynamical core using the horizontal spectral element method and vertical finite difference method: 2-D aspects
نویسندگان
چکیده
The non-hydrostatic (NH) compressible Euler equations for dry atmosphere were solved in a simplified two-dimensional (2-D) slice framework employing a spectral element method (SEM) for the horizontal discretization and a finite difference method (FDM) for the vertical discretization. By using horizontal SEM, which decomposes the physical domain into smaller pieces with a small communication stencil, a high level of scalability can be achieved. By using vertical FDM, an easy method for coupling the dynamics and existing physics packages can be provided. The SEM uses high-order nodal basis functions associated with Lagrange polynomials based on Gauss–Lobatto–Legendre (GLL) quadrature points. The FDM employs a third-order upwind-biased scheme for the vertical flux terms and a centered finite difference scheme for the vertical derivative and integral terms. For temporal integration, a time-split, third-order Runge–Kutta (RK3) integration technique was applied. The Euler equations that were used here are in flux form based on the hydrostatic pressure vertical coordinate. The equations are the same as those used in the Weather Research and Forecasting (WRF) model, but a hybrid sigma–pressure vertical coordinate was implemented in this model. We validated the model by conducting the widely used standard tests: linear hydrostatic mountain wave, tracer advection, and gravity wave over the Schär-type mountain, as well as density current, inertia–gravity wave, and rising thermal bubble. The results from these tests demonstrated that the model using the horizontal SEM and the vertical FDM is accurate and robust provided sufficient diffusion is applied. The results with various horizontal resolutions also showed convergence of second-order accuracy due to the accuracy of the time integration scheme and that of the vertical direction, although high-order basis functions were used in the horizontal. By using the 2-D slice model, we effectively showed that the combined spatial discretization method of the spectral element and finite difference methods in the horizontal and vertical directions, respectively, offers a viable method for development of an NH dynamical core.
منابع مشابه
Verification of a non-hydrostatic dynamical core
The non-hydrostatic (NH) compressible Euler equations of dry atmosphere are solved in a simplified two dimensional (2-D) slice framework employing a spectral element method (SEM) for the horizontal discretization and a finite difference method (FDM) for the vertical discretization. The SEM uses high-order nodal basis functions asso-5 ciated with Lagrange polynomials based on Gauss–Lobatto–Legen...
متن کاملVerification of a non-hydrostatic dynamical core using horizontally spectral element vertically finite difference method: 2D Aspects
متن کامل
CAM-SE: A scalable spectral element dynamical core for the Community Atmosphere Model
The Community Atmosphere Model (CAM) version 5 includes a spectral element dynamical core option from NCAR’s High-Order Method Modeling Environment. It is a continuous Galerkin spectral finite element method designed for fully unstructured quadrilateral meshes. The current configurations in CAM are based on the cubedsphere grid. The main motivation for including a spectral element dynamical cor...
متن کاملBreak fill-dam simulation using GeoStudio (case study: Safarood earth dam, Kerman)
One of the main reason of earth-dam breaks is inner erosion which is caused by washing fine soil particles and finally lead to seepage inside the dam. Seepage will decrease effective tension in the earth-dam's core and arching phenomena is occurred. More settlement of core than to shell is the consequence of arching and occurrence of cracking. In this paper, GeoStudio software based on finite e...
متن کاملA high-order staggered finite-element vertical discretization for non-hydrostatic atmospheric models
Atmospheric modeling systems require economical methods to solve the non-hydrostatic Euler equations. Two major differences between hydrostatic models and a full non-hydrostatic description lies in the vertical velocity tendency and numerical stiffness associated with sound waves. In this work we introduce a new arbitrary-order vertical discretization entitled the staggered nodal finiteelement ...
متن کامل